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A simple method to calculate the static electric polarization of single-walled carbon nanotube 
(SWNT) is obtained within the second-order perturbation approximation. The results are in agree- 
ment with the previous calculation within the random-phase approximation. We also find that the 
low energy electronic band structure of SWNT can be affected by an external electric field perpen- 
dicular to the axis of the tube. The field-induced energy gap shows strong dependence on the electric 
field and the size of the tubes, and a metal-insulator transition is predicted for all (n, n) tubes. Uni- 
versal scaling is found for the gap as function of the electric field and the radius of SWNTs. The 
fact that the external field required to induce a 0.1 eV gap in SWNTs can be reached under the 
currently available experimental conditions indicates a possibility of further applying nanotubes to 
electric signal-controlled nanoscale switching devices. 
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The prospect of nanoscale electronic devices has en- 
gaged great interest. Single-walled carbon nanotubes 
(SWNTs) are significant application as nanoscale de- 
vices due to their extraordinarily small diameter and 
versatile electronic properties It is suggested that 
individual SWNT may act as devices such as field-effect 
transistors (FETs) single-electron-tunncling transis- 
tors ||J|], rectifiers ||J^, or p-n junctions l^j. The 
most exciting expectancy lies in the devices fabricated 
on a single tube In recent years, the interplay be- 
tween mechanical deformation and electrical properties of 
SWNTs have been extensively studied |^-|l2|. Tombler 
et al. ]l^ , p^ used an atomic force microscope tip to ma- 
nipulate a metallic SWNT, leading to a reversible two- 
order magnitude change of conductance, and Lammert et 
al. applied a uniaxial stress to squash SWNTs and 
detect a similar reversible metal-insulator (M-I) transi- 
tion. It is also well known that a magnetic field can also 
change the conductance of carbon nanotubes [|6|-18|. A 
possible electric field-controlled M-I transition are consid- 
ered to be more exciting because of its easy implementa- 
tion in the actual applications. Yet, a question remains: 
Can electric field change the electronic properties of a 
tube? 

In previous studies on electronic transports ||l9|| , the 
potential of a weak longitudinal electric field (bias volt- 
age) in conductors was treated approximately to a slowly 
change variable in the range of the primitive unit cell. 
The electric field makes all the electronic energy and the 
Fermi level have a gradient along the field direction, but 
the energy-band structure is not change. The controlled 
potential, such as the gate voltage without a drop of com- 
ponent in the direction perpendicular to the tube axis in 
the case of FET, is only used to shift the Fermi level or 
changed the carrier concentration In the literature, 
according to our knowledge, there is no report on using 
a transverse electric field to control the longitudinal elec- 



tronic transport of conductors. In (n, n) metallic SWNT, 
the electrons nearby Fermi energy are nonlocal in the cir- 
cumference of the tube since their circumference-Fermi 
wavevector is zero j^], the classic wave-package approx- 
imation in slow-change voltage may be not suitable in 
the presence of the strong transverse electric field. The 
V/nm order electric field is enough to obviously break the 
rotational symmetry about the tube axis, and create new 
interband and intraband coupling, which may change the 
low energy electronic properties of SWNTs, and hence af- 
fect the electronic transport. In the other hand, the field 
is still less than the 2 '-^ 3 orders of the atomic interior 
electric field, can be treated as perturbation. 

In this Letter, we first report the result by using a 
tight-binding (TB) model to calculate the polarizabil- 
ity of SWNTs in the application of an external elec- 
tric field perpendicular to the tube axis. The calcu- 
lated polarizability of SWNT is in agreement with the 
previous results within the random-phase approximation 
(RPA) Then we calculate the low energy electronic 
structure of SWNTs in the electric field. The results show 
obviously valuable effects: (1) The electric field can al- 
ways induce an energy gap in (n, n) metallic SWNTs; (2) 
There is a maximum gap strongly depended on the ra- 
dius of the tubes; (3) Universal scaling is found for the 
gap as a function of the field and the size of the tubes, 
and the numerical results are testified by the second or- 
der perturbation calculations. Our results indicate that 
the magnitude of the electric field required to induce a 
sizable energy gap in metallic SWNTs falls into the range 
of currently available experimental conditions. 

In density function theory, single electronic Kohn- 
Sham Hamiltonian is 



H = T + VKs[p{r)]+Ve{r), 



(1) 



where Vks[p{^)] is the effective potential, self- 
consistently depended on the electronic density p{r), 
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hence the external electric field Eg. Ve(r) is the electro- 
static potential of Eg. Vp{r) = VKs[p{y^)] - Vft:s[po(r)] 
can be acted as the contribution of the polarized charge, 
where po is the unperturbed electronic density. So we 
can rewritten dll) as 



n = T + VKs[poir)] + Vir) 



(2) 



where V{r) — Ve{r) + Vp{r) is the total perturbed poten- 
tial. We have 



Vir) = K(r) + J y^^dv' + A14c(r) 



(3) 



where P6(r) = p(r) — po(r) is the polarized charge density, 
Vxc is the exchange-correlation potential. If V{r) or Ve{r) 
is very small, using the second-order perturbation theory, 
we know AU = ^ J Pb{j^)V{r)dr and 



AU^^J G{r,r')V{r)V{r')drdr' 
AVxcir) = 



(4) 
(5) 



where AU is the change of total energy of electrons, 
G{r, r') can be gained from the unperturbed electronic 
wave functions ^'i(r) and the unperturbed energy levels 
Ei, 



G(r,r') = 



^ ^ E,- E. 

i—occ j—unocc •' 



■[p.,(r)p*,(r')+c.c] (6) 



where py (r) = ^*{r)^j{r), or, 

AC/ 



G(r,r') = 



SV{r)SV{r')' 



(7) 



For an approximation, we treat the V^(r) of as a 
electrostatic potential. If we suppose the polarization in 
SWNT is uniform and due to the induced-surface-charge 
distribution under a constant external field Eg po|| , we 
have E — Ee ~ where p — aoE is the dipole mo- 
ment per length and R being the radius of tubes. By 
selecting a set basis, using slater Xa local exchange- 
correlation approximation, we have numerically solved 
in small tubes with the calculated TB-electron-energy 
change AU. We found the nonhomogenous effect is very 
small . Here we will only report the results within the 
uniform polarization approximation. From — ^ 
we gain the dielectric function e. 



e = Ee/E = 1 



,00 



(8) 



Considering the screened electrons are not confined to 
the surface containing the ions in tubes, after calculat- 
ing the ckQ within RPA, Benedict et al. replaced the R 
of (||) with Reff = R + SR for calculating e. Fitting 
the result of Cgo, they gave SR — 1.2 A for some small 



radius tubes. If Directly starting from the exactly un- 
perturbed wave functions and energy levels of tubes, we 
can gain the real distribution of polarized charges. By 
numerically solving (||), (||) and (0), we found that the 
exchange-correlation potential maybe decrease partly the 
depolarized field which is completely attributed to Reff 
by Benedict et al. However, the order of magnitude of 
e which we need can be identically obtained from one of 
both (|) and (|). 

The nearest-neighbor TB Hamiltonian has been used 
successfully for calculating the electronic structure of 
graphite sheet and nanotubes 0. Our previous 
works show the model can well describe the elec- 
tronic perporties and total energy. The Hamiltonian is 
dependent on the matrix element of potential V{r), 

<4(r-RO|7^|0,(r-Rfc) >^ Hi' {l,k) + V,,{l,k), (9) 

where |0i(r — Hi) > is the ith (2s or 2p) orbit of the R; 
carbon atom, TCq {I, k) are the Slater-Koster TB param- 
eters in the absence of electric field |22j . Expanding the 
V{y) near the lattice vector R;, 



F(r + RO«l^(Ri)-E(RO-r, 



(10) 



where we suppose the E is uniform or slowly change in 
the size of atom. Only considering the nearest-neighbor 
contribution, we have 

V,j{l,l) « %F(RO - E(RO • P,,(0), (11) 
F,,(/, k) « s.,l-(R?J - E(RJ=,.) . P,,(d), (12) 

where Pij(O) and (d) are the dipole matrix elements, 
and Sij are overlap integrals which are only slightly affect 
the electronic energy bands, d = R; — R^, and R^j, = 
(R; -\-TIk)/2 is the position of mass centre of / atom and 
k atom. Benedict et al. showed that the P^ can be 
neglected except the Psp(O) — RgpX « 0.5x A. From 
the Hamiltonian, we can calculate the total electronic 
energy in presence of any total electric field E(r). Since 
the polarizability of a electrons is smaller than that of tt, 
we only calculate the contribution of tt. The results are 
shown at Fig.|[ The values of ao of a few small radius 
tubes are in agreement with the previous results Q , but 
aQ « 2.4_R2-4 jg faster increase as the radius increase than 
the result of Benedict et al. (ao ~ 2.6i?2). Fig.||(b) shows 
the obtained value of e of some (n, n) tubes using the R 
and Reff, respectively. The e is about order of 5 ~ 15 
for a few small tubes, slowly increases to about 20 for 
(60, 60) tube. Since the polarization is nearly constant 
in SWNT, we will only calculate the effects of a uniform 
total electric field in this letter. Fig.|| shows the electronic 
energy bands of a (10, 10) tube. When Vq = 0.5 V, where 
Vq = ER, a sizable gap about 0.3 eV is found at Kq, the 
Fermi wavevector in zero field. As E increases, the gap 
increases, and when Vq = 1.5 V, the bands structure 
is obviously deformed. It is surprising to find that the 
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gap decreases as E increases further. When Vq — 3.0 V, 
the zero gap is found, but the Fermi point dramaticaUy 
moves from Kq. 

To probe the above effect in general, we performed 
the computation for a series of (n, n) tubes. Fig.^ shows 
the gap as a function of the apphed field Vq in (n, n) 
tubes, where n is from 5 to 15. From the figure, we 
find the determined effect: The electric field can always 
induce a gap in (n, n) tubes, and the size of the gap 
strongly depends on the amplitude of the transverse field 
and the tube parameter n. For any (n, n) tubes, the 
gap first increases with increasing field, and reaches a 
maximum value Egm at the Vom, then drops again. Both 
the maximum gap Egm and the corresponding Vom are 
approximately proportional to 1/n, and hence inversely 
proportional to the radius of tubes, i.e., 

Egm ~ 6.89 eV/n, e„ w 12.09 eV/n. (13) 

The finding shown in ( |l3| ) that electric field effects fastly 
decrease as the size of conductor increases might be the 
reason why people have not yet recognized the effect in 
previous studies. The field dependence of the gap is quite 
similar for tubes with various cross-sectional radii, which 
invokes us to scale both Eg and Vq up n times their orig- 
inal values. The obtained results are shown in Fig||. 
From it we do find the scaled gap to be a universal func- 
tion of the scaled electric field for all {n, n) tubes. In 
the low field range, for all calculated eleven (n, n) tubes 
there exists a simple relation: nEg = A(nVb)^, where 
A is a constant, about 0.07 (eV)~^. In the higher field 
range, except for a few small-radius tubes such as (5, 5) 
and (6, 6) tubes, the universal scaling law still holds. 

To understand the above scaling relation, we use per- 
turbation theory to calculate the field-induced gap in low 
field limit. The first-order perturbation approximation 
only causes shift in the Fermi level, showing no contri- 
bution to the gap change. Calculating up to the second 
order perturbation at Kq point, we obtained the follow- 
ing analytic result 

E,^^^nEV. ,14) 

where h (=3.033 eV) is the hopping parameter in the ab- 
sence of the electric field |2^,Q . The contribution of the 
overlap integral s, which is very small, is neglected. Ob- 
viously, the second-order perturbation calculation gives 
almost the same scaling relation as the numerical results 
in the low field, though the obtained A « 0.09 (eV)^^ is 
slightly larger than the numerical result 0.07 (eV)^^. In 
the high field range, since the Fermi wavevector is moved 
from Kq, the perturbation theory becomes not suitable. 
However, the low field range may be more compatible 
with the practical application. In order to open a 0.1 eV 
gap in the energy bands of {n, n) tubes, n must be smaller 
than no = 6.89/£'g ~ 68, and the required electric field 
is. 



\E\-^J^n-K (15) 
Stqc V A 

where vq (=1.42 A) is the bond length of carbon atoms 
in the SWNT. Therefore, for example for a (10, 10) tube, 
the required field is about 5 x 10^ V/m, and for a (60, 60) 
tube, it is about 3x 10^ V/m. Considering the polarizing, 
the needed external field is about 0.5 V/A and 0.06 V/A, 
respectively. Our results are encouraged by an important 
fact that the perturbation approximation is suitable only 
required the small total field E, need not the small exter- 
nal field Eg = eE. Even though e is larger, our conclusion 
is still correct, but needing stronger external field. The 
magnitude of the required external field for inducing a 
sizable gap in tubes with larger radius can be reached by 
the currently available experimental conditions, we wish 
the above prediction can be checked in near future. 

In summary, we have proposed an electric field-induced 
M-I transition in (n, n) SWNTs for the first time. The 
results support the argument that SWNTs can be applied 
as nanoscale electric signal-controlled switching devices. 
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sions on the results. The numerical calculations were 
performed partly at ITP-Net and partly at the State Key 
Lab. of Scientific and Engineering Computing. 
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Fig. 1 X. Zhou et al 

FIG. 1. A uniform external electric field perpendicular to 
the axis of SWNT. 
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Fig.2 X. Zhou et al 



FIG. 2. The polarization of {n,n) tubes: (a) ao ~ vP , 
where 7 w 2.4 (see text); (b) the dielectric function of tubes, 
circles and squares are the calculated results using R and 
Reff, respectively (see text). 
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Fig. 3 X. Zhou et al 

FIG. 3. The energy bands of (10, 10) tube in the vicinity 
of the Fermi level under the application of a transverse elec- 
tric field of certain magnitude. Vb = OV is the result in the 
absence of the electric field. 
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Fig. 4 X. Zhou et al 

FIG. 4. Field-induced gap of (n, n) tubes versus the field 
Vb = ER. From the top to the bottom, the tube parameter n 
increases from 5 to 15. The clubs denote the position of the 
maximum gap point. In the inset, both the maximum gap 
Egm and its corresponding electric field Vbm are found to be 
proportional to 1/n The lines are fitting results. 
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Fig. 5 X. Zhou et al 

FIG. 5. Universal scaling are found in the gap as a function 
of Vb in different (n, n) tubes. In the low field, the data of 
all tubes are very much consistent with the scaling relation 
nEg = \{nER)^, as expected by the second perturbation 
theory. The line is the fitting result. Except for (5, 5) and 
(6, 6) tubes, the universal scaling is satisfied well up to high 
field region. 
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